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ABSTRACT 
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. , Context. The rate of star formation both in the Galaxy and in external galaxies should be related to the physical properties of the 

\^ . molecular clouds from which stars form. This is expected for the starbursts found both in irregular galaxies and in some mergers. The 

dwarf galaxy Henize 2-10 is particularly interesting in this context as it shows a number of newly formed Super Star Clusters (SSCs) 
associated with a very rich molecular environment. 

Aims. We present a high angular resolution study of the molecular gas associated with the SSCs with the aim of deriving the physical 
^ . properties of the parent molecular clouds. The final goal is to test the expectation that the formation of SSCs requires exceptionally 

( ' dense and massive clouds. 

lypj ' Methods. We have used the Submillimeter Array with an angular resolution of l'!9 x I'.'S to map the J=2-l transition of CO in 

, Henize 2-10. Supplementary measurements of HCN(y=l-0), '^^00(7=2-1) and millimeter continuum were obtained with the APEX, 

IRAM 30m and SEST single dish telescopes. 

Results. Our single dish observations confirm the association of the newly formed SSCs in Henize 2-10 with dense molecular gas. 
Our interferometric observations resolve the CO(2-l) emission in several giant molecular clouds. Overall the molecular gas accounts 
^ I for approximately half of the mass in the central regions of Henize 2-10. Although we find indications that the molecular clouds 

\l . associated with the formation of SSCs in Henize 2-10 are massive and dense, the tracer we used (CO) and the linear resolution of our 

observations (60x80 pc) are still not adequate to test the expectation that exceptionally dense and massive cores are required for SSCs 
formation. 

o 

1. Introduction Orion cloud neighbouring the ONC. Moreover, both gas pres- 

sures and star formation rates are higher in the inner regions 
Galactic and extragalactic star formation mostly occurs in of the Milky Way, where most of the young massive star clus- 
"bursts" or clusters and the majority of stars is thought to form jg^s are located, and in external galaxies (e.g. Kennicutt 1998). 
in large clusters. However, cluster masses and densities vai-y ^ summary of the characteristics of both clouds and clusters is 
d gready even within the Milky Way (see Lada & Lada 2003) given by Tan (2008) 
where the well studied Orion nebular cluster (ONC) has been 

estimated to have a mass of order 1800 Mq (Hillenbrand & Of particular interest are the youngest extragalactic embed- 
Hartmann 1998) while the most massive galactic clusters (e.g. ded super star clusters or SSCs which have been identified by 
Brandner et al. 2008; Dowell et al. 2008) have masses of a the free-free radio emission from their associated HII regions 
few times 10"* Mq. These are dwarfed however by young extra- (e.g. Kobulnicky & Johnson 1999; Tarchi et al. 2000; Johnson 
galactic clusters and in particular by the starbursts or super star & Kobulnicky 2003). These are similar in many ways to local 
clusters (SSCs) seen in merging systems (e.g. Whitmore 2002) ulti-a-compact HII regions but larger and ai-e typically powered 
and in some iiTegular galaxies (e.g. Elmegreen 2002; Johnson & by the equivalent of several thousand O stars (e.g. Kobulnicky 
Kobulnicky 2003) where masses can reach lO*" Mq. It has been & Johnson 1999). It is estimated that their stellar masses exceed 
speculated that in this case one may be observing young globular ~ 10^ Mq, their radii are of order of a few parsecs, their elec- 
clusters. tron densities are of order of lO'* cm"^ and their ages are below 

One presumes that the diversity between local and distant 1 Myr (Kobulnicky & Johnson 1999; Vacca et al. 2002; Johnson 
star formation reflects variation within the molecular clouds & Kobulnicky 2003). Thus their mass is of the same order as 
from which these clusters form and, in particular, variations in galactic GMCs but their dimensions are smaller by an order of 
the gas pressure. Indeed, even locally, one observes large dif- magnitude and, correspondingly their densities and pressures are 
ferences in the characteristics of nearby clouds with moderate much larger. For an assumed star forming efficiency of around 
pressures and densities in Taurus cores, which mostly form iso- 10%, we might expect to see molecular forerunners of the SSCs 
lated stars, and higher pressures in the molecular gas in the with masses of 10^ Mq and sizes of a few parsecs. One might 
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hope to detect such precluster "cores" either in molecular line 
emission or in continuum dust emission. 

The dwarf irregular galaxy Henize 2-10 at a distance of 9 
Mpc (i/o=75 km s"' Mpc'^ JVacca & Conti 19921 Kobulnicky & 
Johnson 1999) is a good example of this phenomenon. It con- 
tains several times 10^ Mq of molecular gas as well as multiple 
HII regions and young clusters observed at radio, optical, and 
infrared wavelengths (Kobulnicky & Johnson 1999; Johnson et 
al. 2000; Vacca et al. 2002; Cabanac et al. 2005). The nebulosity 
seen in HST images does not coincide well with the radio knots 
suggesting large amounts of obscuration. It is clear therefore that 
higher angular resolution observations of the molecular emis- 
sion are warranted in order to understand the characteristics of 
the clouds from which the SSCs form. Ultimately of course, one 
would like to image molecular line emission with a resolution 
comparable both to the dimensions of the SSCs and the optical 
images (0.1 arcsec resolution). This is currently not possible but 
as a first step in this direction, we present here results obtained 
with the Submillimeter array (SMA) in the CO J=2-l transition 
(and in the mm continuum) with a resolution of roughly I'.'S 
(80 parsec at the distance of Henize 2-10). This allows us to see 
the distribution of structures comparable to galactic GMCs and 
to compare with images at other wavelengths. It also allows us 
an improved view of the molecular cloud dynamics than avail- 
able in earlier work (Kobulnicky et al. 1995). Finally, we have 
obtained supplementary measurements of the millimeter con- 
tinuum and the HCN(y=l-0) and '^CO(7=2-l) lines with the 
SEST, IRAM 30m and APEX single dish telescopes. 

This paper is structured as follows. In section |2] we present 
our observations. In section [3] we present our results and com- 
pare with previous CO observations of Henize 2-10. In section 
m we describe our analysis and derive the physical parameters 
of the detected clouds. In section |5j we discuss the results and 
finally, in section|6l we present our summary and conclusions. 



2. Observations and Data Reduction 

2.1. Single-dish observations 
2.1.1. I RAM-30m Telescope 

Henize 2-10 was observed in the ^^CO(J^2-l) (220.399 GHz) 
and HCN(7=l-0) (88.63 GHz) Unes on December 2007 with 
the 30m Institut de Radioastronomie Millimetrique (IRAIVQ) 
telescope on Pico Veleta. A single position (08*^36'" 15^.20, - 
26°24'34'.'0 [J2000]) was observed. The spectra have a typical 
beam size of -11" in the "COiJ^2-l) line and -28" in the 
HCN(7=l-0) Une. Calibration, pointing and focus checks were 
performed regularly. 

The data processing was done with the CLAS^I softwares. 
The data were converted into flux density units, using S/T^b - 
4.95 Jy/K (Rohlfs & Wilson 1996). The ^^CO(J^2-l) and 
HCN(7=l-0) spectra were smoothed to a velocity resolution of 
5 km s~^ and 10 km s respectively, and the final rms noise was 
40 mJy and 5 mJy per channel, respectively. 



' This publication is partially based on observations carried out with 
the IRAM 30m telescope. IRAM is supported by INSU/CNRS (France), 
MPG (Germany) and IGN (Spain). 

^ See http : //w vy w. iram.fr/IR AMFR/GILDAS] for more information 
about the GILDAS softwares. 



2.1.2. APEX Telescope 

We observed Henize 2-10 in the ^^CO(J=2-l) (220.399 GHz) 
with the Atacama Pathfinder Experiment telescope (APEXlJ, 
Gusten et al. 2006) on August 28th and 31st 2008 and on 
December 17th 2008, in the same position used for the IRAM- 
30m observations. The typical FWHM beam size for the 12m 
antenna at the '^CO(y=2-l) frequency is 28". The line was 
measured in upper sideband (USB) using the APEX-1 facil- 
ity receiver in combination with the Fast Fourier Transform 
Spectrometer. 

The data were reduced using the CLASS software. A first 
order polynomial was fitted to the line-free channels and sub- 
tracted off' the baseline. We have converted our line intensities 
into a flux density scale, assuming the aperture efficiency obser- 
vationally determined by Gusten et al. (2006): 1 K of antenna 
temperature corresponds to 39 Jy. To increase the signal to noise 
ratio, the spectrum was smoothed to a velocity resolution of 10.7 
km s The noise level in our final spectrum is of the order of 
70 mJy per velocity channel. 

2.1.3. SEST Telescope 

Henize 2-10 was mapped in the 1.2 mm continuum with the 
bolometer array SIMBA at the SEST on the nights 7-8 and 8- 
9 June 2002. The beam size of the SEST at this wavelength 
is 21". Background subtraction was obtained by fast scanning 
(80"/sec) of the source through the array field of view in AZ and 
ALT, the map size was 480" x 240". The atmophere opacity was 
T ~ 0.1 during the observations. We reduced the data with the 
MOPSI (Zylka 1998) package. The total integration time was 4 
hours and 50 minutes. The flux calibration was obtained by com- 
parison with planet Uranus. We detected an un-resolved source 
whose flux is 52+5 mJy. 

2.2. Interferometric SMA observations 

Submillimeter interferometric observations of Henize 2-10 in 
the CO(7=2-l) (230.538 GHz) line were carried out on February 
and March 2008 with the Submillimeter ArrajQ (Ho et al. 2004). 
We obtained data with the array in the extended and compact 
configurations, half track each, and with six antennas and seven 
antennas respectively. The phase centre of the observations was 
the same position used for the IRAM-30m observations. The 
range of baselines of our dataset spans from the shadowing limit 
to ~ 220 m. The observations cycled between Henize 2-10 and 
two gain calibrators (0826-225 and 0730-116), with 15 minutes 
on the target and 7 minutes on the calibrators. Passband calibra- 
tion was ensured observing 3C273 and 3C279. 0730-116 was 
used to set the absolute flux scale. The flux of this calibrator at 
the time of observations was determined using values from mon- 
itoring observations of the SMA standards (moons and planets). 
The SMA receivers operate in a double-side band mode, with 
the upper and lower side bands separated by 8 GHz. We tuned 



^ This publication is partially based on data acquired with 
the Atacama Pathfinder Experiment (APEX, ESO project 081.F- 
9813). APEX is a collaboration between the Max-Planck-Institut fur 
Radioastronomie, the European Southern Observatory, and the Onsala 
Space Observatory. 

The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 
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to observe the CO(2-l) in the upper sideband and the '"'CO(2-l) 
and C'^0(2-l) in the lower sideband. 

The initial flagging and calibration of the raw visibility data 
were done using the IDL-based MIR packagfl The MIRIAD 
package was used for imaging and deconvolution of the cali- 
brated visibilities and for the data analysis. The continuum was 
constructed from the line free channels in the visibilitiy do- 
main. Continuum and line maps from the combined datasets 
were made with natural weighting, providing a final synthesized 
beam FWHM of r.'9xr.'3 (at a position angle of 5 deg) for the 
line map, which corresponds to a physical size of 80x60 pc for 
Henize 2-10, and of r.'9x r.'4 (at a position angle of 6 deg) for the 
continuum map. The rms noise is 19 mJy beam"' per 5 km s"' 
channel in the CO(7=2-l) line data and 1 mJy beam ' in the 
continuum. 



3. Results 

3.1. Millimeter Continuum 

No millimeter continuum detection was obtained at the high an- 
gular resolution of the SMA. We obtain no clear detection above 
a 5 cr upper limit of 5 mJy/beam, which converts into a total 
gas-Hdust mass between 5 10^ Mq and 1.6 10^ Mq, assuming a 
range of temperatures between 20 K and 50 K and a dust opacity 
of 0.005 cm- g"' (see discussion in section l431 l. 

Our SIMBA single dish 1.2 mm flux is consistent with the 
MAMBO 1.3 mm flux of - (56 ± 14) mJy reported by Galliano et 
al. (2005), while, integrating the SMA continuum map over the 
30m beam, we recover 17 mJy with the interferometer, which 
is the 30% of the single dish flux. The most likely explanation, 
given also the extent (~15") of the 450 and 850 jim emission 
detected by Galliano et al. (2005), is that our interferometric ob- 
servations filter out the most extended emission. To estimate this 
effect we used the approach outlined by Wilner & Welch (1994). 
Given a minimum baseline length, their formalism allows us to 
estimate the fraction of flux seen by the interferometer for a 
Gaussian source with a given FWHM. For our observations, the 
minimum projected baseline of 7 kA, where A is the wavelength 
of our observations, allows us to recover 100% of the total flux 
for sources comparable to the beamsize, but only ~40% of the 
total flux density for a source of FWHM 15" (650 pc at a dis- 
tance of 9 Mpc). 



3.2. Molecular Gas Morphology 
3.2.1. Single dish observations 

In Fig. [T] we show the observed IRAM-30m spectra of the 
'^CO(7=2-l) and HCN(7=l-0) transitions. Both lines are de- 
tected and they are approximately 100 km s"' wide, consis- 
tent with previous observations of CO and isotopologues by 
Kobulnicky et al. (1995). 

The line parameters have been estimated fitting a Gaussian 
profile to each spectrum. Table [1] summarizes the beam size at 
the observed frequencies, the peak and integrated flux densties 
and the line width of the observed emission lines. 

The HCN(l-O) total intensity that we derive from our 
spectrum (<2 Jykms"') is significantly lower than the value 
found by Imanishi et al. (2007) in their NMA integrated map 
(~6 Jykms '). This is surprising as the single dish telescope 



800 



900 




0.02 



-0.02 



0.02 



800 900 

LSR velocity (km/s) 
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Fig. 2. Comparison between the '^CO(7=2-l) spectra from 
IRAM-30m (grey shaded), SMA (solid line, integrated over the 
IRAM-30m beam area), JCMT (dotted line; Baas et al. 1994) 
and APEX (dashed Une). 



should have detected the same or more emission than the inter- 
ferometer Both the IRAM-30m and the NMA detections do not 
have a very high signal to noise ratio and higher sensitivity ob- 
servations will be needed in the future to confirm the flux of this 
source. We cannot exclude the possibility of a problem with the 
calibration of the IRAM spectrum. 

In Fig. |2]we show the observed APEX '^CO(y=2-l) spec- 
trum (dashed line). The parameters of the emission line from a 
Gaussian fit to the spectrum are presented in Table[T] In the same 
figure we show also the IRAM-30m spectrum (grey shaded) and 
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Table 1. IRAM-30m and APEX '^CO(y=2-l) and HCN(7=l-0) emission lines parameters. 
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Fig. 3. Comparison between the NRAO-12m CO(7=2-l) spec- 
trum (dashed line) from Kobulnicky et al. (1995), our SMA 
CO(y=2-l) spectrum (solid line), integrated over the whole re- 
gion of the emission, and the JCMT CO(7=2-l) spectrum (dot- 
ted Une) from Baas et al. (1994). 

the JCMT spectrum from Baas et al. (1994, dotted line). Figure 
|2] shows that the IRAM-30m '^CO(/=2-l) spectrum is fainter 
than both the JCMT and the APEX spectra. This could be due to 
either a problem in the calibration or pointing of the IRAM-30m 
data, because of the large airmass of this southern object, or to 
the fact that the '^CO(7=2-l) emission is more extended than 
the IRAM-30m beam (11" at this frequency). Therefore part of 
the emission which is seen by APEX and JCMT, with beam 
sizes of 28" and 21" respectively, is outside the IRAM beam. 
This explanation is consistent with the SCUBA continuum maps 
of GalHano et al. (2005), the non-detection of the '^CO(7=2-l) 
emission with the SMA (see below) and the detection of an ex- 
tended component in the CO(7=3-2) emission by Vanzi et al. 
(2009). 

3.2.2. Comparison between single disli and interferometer 

Figure [3] shows the NRAO CO(y=2-l) spectrum (dashed line) 
from Kobulnicky et al. (1995) overlaid on our SMA integrated 
spectrum (solid line) and the JCMT CO(7=2-l) spectrum from 
Baas et al. (1994, dotted line). The comparison shows that we 
recover 44% of the single dish flux density, while the remain- 



ing part, which most probably comes from extended emission, is 
filtered out by the interferometer. 

Kobulnicky et al. (1995) found an elongated cloud of molec- 
ular gas with a tail towards the south-east. This tail, revealed 
in the OVRO CO(7=1-0) map in a velocity interval between 
805 and 835 km s"', is not seen in our SMA CO(7=2-l) data 
(with the exception of two 3 cr level peaks), most likely because 
the higher resolution SMA observations filter out this extended 
emission. 

As mentioned before, the correlator frequency coverage and 
configuration were such that also '^CO(7=2-l) and C'^0(7=2- 
1) transitions were observed. From our IRAM-30m observations 
of '''CO(/=2-l) we expect a peak flux density of about 0.2 
Jy on a beam size of around 11". The spectrum of our SMA 
'^CO(7=2-l) observations integrated over the IRAM beam is 
shown in Fig.|2] overlaid on the IRAM-30m '^^CO(y=2-l) spec- 
trum, the JCMT '^CO(y=2-l) spectrum from Baas et al. (1994) 
and the APEX '^CO(y=2-l) spectrum. The rms of the SMA 
spectrum (roughly 0. 1 Jy) is about half of the peak flux expected 
from the single dish IRAM observations and does not permit 
a credible detection. The C"*0(7=2-l) emission, as one might 
expect, was detected neither in the APEX spectrum nor by the 
SMA. 

The non-detection of the '^CO(2-l) line and the large frac- 
tion of the CO(2-l) flux that is not recovered by our SMA obser- 
vations are consistent with the presence of extended molecular 
emission to which the SMA is not sensitive, as discussed in the 
previous paragraph, and with the fact that the SMA is only re- 
covering a fraction of the millimeter continuum emission. 

3.2.3. IVlorpliology of tine molecular gas emission 

Our SMA spectral Une data resolve the OVRO CO(7=1-0) emis- 
sion from Kobulnicky et al. (1995) into several compact sources. 
The velocity structure of the CO(7=2-l) emission is given in 
Fig. m which shows the channel map of the CO(7=2-l) emis- 
sion. The map reveals the complex velocity structure of the CO 
emission, with most of the emission centred on an LSR velocity 
around 850 km s"'. Several velocity components can be distin- 
guished in the emission and some of them are comparable with 
the beam size. The CO components are well separated in both 
position and velocity and spread over a velocity interval of about 
100 km s as already revealed from the single dish observations 
of '^CO(y=2-l) and HCN(y=l-0) (see Fig.Qand Fig.©. 

Figure |5] shows the SMA CO(7=2-l) velocity-integrated in- 
tensity map in solid contours overlaid on the VLA 3.6 cm con- 
tinuum (Johnson & Kobulnicky 2003). Part of the centimeter 
continuum emission, associated with the super star clusters, is 
coincident with the most compact molecular gas sources in the 
CO(y=2-l) emission line. In particular, the centimeter sources 
1, 2 and 5 (see Johnson & Kobulnicky 2003) correspond respec- 
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Fig. 4. Channel map of the CO(7=2-l) line. The contour levels are from 5 cr {o--l9 mJy/beam) in steps of 3 cr. The black numbered 
crosses indicate the peak position (in the channel corresponding to the peak velocity) of the emission of the 14 molecular clouds 
identified by the Clumpfind algorithm (see ^4.1l and TableO. 




Table 2. Integrated SMA CO(7=2-l) and "CO(7=2-l) fines pa- 
rameters. 



R.A. 

Fig. 5. SMA integrated CO(7=2-l) map (contours), overlaid on 
the VLA 3.6 cm continuum image from Johnson & Kobulnicky 
(2003). The contour levels are from 5 cr in steps of 2 cr (cr=0.8 
Jy beam ' km s '). The white line and tfie cross indicate respec- 
tively tfie orientation of tfie position-velocity cut, sfiown in Fig|6] 
(see ^5.1l l. and tfie central position of tfie map. 
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tively to our sources 12, 13 and 6 (see Fig.|4]and TableO. Some 
of tfie CO emission is not associated witfi tfie SSCs. 

Tfie strongest 3.6 cm source, wfiicfi does not correspond to 
any significant amount of molecular emission, is likely in an 
older evolutionary stage, fiaving already dispersed most of tfie 
surrounding molecular material from wfiicfi it formed. 

We performed Gaussian fits to tfie spectra of tfie lines ob- 
served witfi SMA, integrated over tfie wfiole region of tfie emis- 
sion, to estimate tfie peak and integrated flux density and tfie 
line widtfi of tfie observed emission lines. Table |2] summarizes 
tfie fitting parameters. Tfie values of tfie peak and integrated flux 
densities of tfie '^CO(7=2-l) emission line, observed witfi SMA, 
are upper limits because tfie line is not detected. It is tfius pos- 
sible to derive a lower limit to tfie flux density ratio between 
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the CO(7=2-l) and '^CO(7=2-l) emission lines. We estimate a 
CO(7=2-l)/i3CO(y=2-l) ratio > 14.9. 

4. The structure of the CO(7=2-1) emission 

4.1. Identification of the clouds: Clumpfind 

We applied the Clumpfind algorithm dWilliams et al. 1994] ) to 
the SMA CO(7=2-l) data cube to decompose the emission 
into discrete structures and identify the individual molecular 
clouds. Clumpfind works by making a contour map through 
the data cube to locate the peak flux densities and then steps 
down the contour levels to determine where the clumps are lo- 
cated. Clumpfind allows the user to adjust the contour spacing 
and the lower contour level. The advantage of this algorithm is 
that we make no assumptions about the cloud shape. The spe- 
cific intensity contour used as emission threshold, which defines 
the borders (boundaries) of the individual identified clouds, is 
0.1 Jy/beam, corresponding to 5 cr level. The step level used 
to disantangle the emission and identify isolated peaks is 60 
mJy/beam (3 cr). To be sure the algorithm is not introducing arte- 
facts, we decided to limit the analysis to the clouds with peak 
flux density larger than 15 cr. In this way, we identified 14 sepa- 
rate CO clouds. 

Is is important to realize that the Clumpfind algorithm does 
not take into account effects of self-absorption, which manifest 
in a double peak line profile: it would then try to separate the two 
peaks of a single self-absorbed clump spectrum in two distinct 
clouds with same spatial position and different velocity. With our 
data we do not see a clear evidence that this is effecting a signif- 
icant fraction of the clumps we identify, as most of the clumps 
are also spatially distinct. 

To emphasize the complex substructure of the emission we 
marked with black crosses in Fig. |4] (the channel map of the 
emission) the position of the clouds identified by the Clumpfind 
algorithm, in the channel relative to the peak of each cloud emis- 
sion. The 14 components are well visible in the plot and are la- 
beled by their respective number. 

The parameters for our final cloud decomposition of the data 
set, given by Clumpfind, are presented in Table [3] The angular 
sizes have been computed from the 5 cr contour levels and de- 
convolved with a simple Gaussian deconvolution for the beam 
size, assuming source and beam to be Gaussian. The velocity 
widths have been deconvolved as well for the velocity resolution. 
Table[3]gives also the derived masses of the identified clouds (see 
next section for the details). 

4.2. Derivation of Physical Parameters of the Clouds 

Several methods can be used to obtain molecular mass esti- 
mates. The molecular hydrogen column density, N(H2), can be 
estimated from the CO flux, ^co, using a CO-to-H2 conversion 
factor, Xco (Xco=N(H2)/Ico), which is empirically determined 
from Galactic molecular clouds. It is a matter of debate what is 
the best choice for Xco- Dwarf galaxies, possibly due to their 
low metallicities, seem to have a higher standard conversion fac- 
tor than large metal-rich spirals (Maloney & Black 1988; Verter 
& Hodge 1995; Wilson 1995; Arimoto et al. 1996). However, 
Henize 2-10 is a dwarf galaxy with roughly solar metallicity 
(12H-log(0/H)=8.93; Vacca & Conti 1992). Moreover, BHtz et al. 
(2007), studying extragalactic Giant Molecular Clouds (GMCs), 
conclude that there is no clear trend of Xco with metallicity. 
We thus estimate the molecular mass (M^oi) of the clouds using 
the Galactic value for the CO-to-H2 conversion factor, Xcogai - 



(3+l)xl02° cm-2 (Kkms-i)-' (Strong etal. 1988; Sanders et 
al. 1987) which also includes a factor of 1.36 to account for the 
helium mass contribution. Following Wilson & Scoville (1990) 
we use: 

M„„, = 1.61xlO^(iH^)'<p,|2 (1) 

where Sco is the CO(y=2-l) flux in Jy km s"', c/Mpc is the dis- 
tance and /?2i is the CO(7=2-l)/CO(y=l-0) line ratio, which we 
assumed equal to 0.89 (IBraine et al. 1 993). The latter value is 
also consistent with the range of ratios derived for this galaxy 
(Meier et al. 2001; Baas et al. 1994) and for dwarf irregular 
galaxies (Petitpas & Wilson 1998). 

The masses of the individual clouds can also be calculated 
from the measured sizes and velocity dispersions using the 
virial theorem. The underlying assumption is that the internal 
kinetic energy and the gravitational energy are in equilibrium 
(Solomo n et al. 19871 1. Assuming the source to be spherical and 
homogeneous and neglecting contributions from magnetic field 
and surface pressure, the virial cloud mass is given by: 

MviR = 0.509 d(kpc) 0s(arcsec) AV^^2(km s"') Mg (2) 

(MacLaren et al. 1988), where AV1/2 is the measured full-width 
velocity at half-maximum intensity in km s^' and ©s is the cloud 
angular diameter in arcsec. 

The masses derived for each cloud by means of the two dif- 
ferent methods described above are presented in Table [3] The 
values show that the two diff'erent methods give very similar 
results (within a factor typically of 2). In particular, the virial 
masses are slightly larger than the masses computed using the 
conversion factor The two methods used here to determine the 
masses of the clouds are in reasonable agreement and no clear 
evidence for a systematic diff'erence between the two determina- 
tions can be claimed at this point. There has been much discus- 
sion of how well CO intensity can be converted to total molecu- 
lar gas mass in external galaxies using Galactic conversion fac- 
tors. However, in the present study we are resolving individual 
molecular clouds rather than averaging over the entire distribu- 
tion of clouds. The fact that the masses derived by the two dif- 
ferent methods agree within a factor of 2, in spite of the diff'er- 
ent assumptions made for the two determinations, suggests that 
there may not be a significant problem in this case. 

We also computed the surface density, 2=M/(7rR^), of each 
identified cloud, using both the virial masses and the masses de- 
rived with the conversion factor The values are given in Table |3] 
(Col. 11 and 12). 

4.3. Continuum emission from the molecular clouds 

It may appear surprising that we do not detect 1.3mm contin- 
uum emission from the clumps seen in the CO emission or from 
the compact HII regions seen with the VLA (see Fig. |4]and|5]l. 
As far as the HII regions are concerned, we note that extrapo- 
lation to 230 GHz of the Johnson & Kobulnicky (2003) flux of 
the strongest centimeter continuum source, assuming optically 
thin free-free emission, yields expected flux densities of 2 mJy 
or roughly 2 cr at our sensitivity. Thus it is not surprising that 
our observations do not detect the free-free emission from these 
source, especially if partially resolved. 
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Vpeak 
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Sc0(J=2-l) 


MviR 




^VIR 






R A rT9nnni 

Iv.rt. [J Z.UUUJ 


DppI rT?nnfii 


[km/s] 


[Jy/beam] 


[arcsec] 


[km/s] 


[Jy km/s] 


[xIO^Mg] 


[xlO'^Mol 


[g/cm-] 


[g/cm^] 


(1) 




CW 
\-') 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


S-^A- 1 '> 
o. jO. iJ.JJ 




840 


0.48 


1.3 


18 


10.7 


3.6 


3.9 


0.08 


0.09 


2 


S-^A- 1 
o.J\). 1 J.J^ 




835 


0.47 


1.2 


25 


13.4 


6.9 


4.9 


0.16 


0.11 


3 


S'^A- 1 

o. jO. i J.jU 


-ZD.ZM-. J 1 .Z 


835 


0.46 


1.2 


18 


7.9 


3.3 


2.9 


0.09 


0.08 


4 




-ZD.ZM-. J / .U 


850 


0.46 


0.8 


19 


7.5 


2.9 


2.7 


0.14 


0.14 


5 


8:36:15.24 


-26:24:32.4 


850 


0.43 


1.2 


25 


12.7 


7.1 


4.6 


0.17 


0.11 


6 


8:36:15.31 


-26:24:35.0 


855 


0.42 


1.0 


20 


8.3 


3.9 


3.0 


0.13 


0.10 


7 


8:36:15.14 


-26:24:37.6 


845 


0.42 


1.2 


13 


8.1 


1.9 


3.0 


0.05 


0.07 


8 


8:36:15.43 


-26:24:36.2 


825 


0.33 


0.8 


20 


5.0 


3.0 


1.8 


0.16 


0.09 


9 


8:36:15.38 


-26:24:36.6 


845 


0.32 


1.2 


26 


8.1 


7.4 


3.0 


0.19 


0.07 


10 


8:36:15.10 


-26:24:34.4 


865 


0.31 


1.1 


20 


4.1 


4.0 


1.5 


0.12 


0.05 


11 


8:36:15.54 


-26:24:37.0 


820 


0.29 


1.0 


15 


5.0 


2.1 


1.8 


0.08 


0.06 


12 


8:36:15.03 


-26:24:33.4 


895 


0.29 


1.0 


19 


5.2 


3.1 


1.9 


0.11 


0.07 


13 


8:36:15.09 


-26:24:34.0 


875 


0.28 


0.9 


19 


5.0 


3.1 


1.8 


0.12 


0.07 


14 


8:36:15.19 


-26:24:37.6 


825 


0.28 


1.2 


24 


5.3 


6.5 


1.9 


0.16 


0.05 



NOTES. - Col.(l): Number of the cloud. Col.(2)-(3): Peak position of the emission. Col. (4): Peak velocity of the emission. Col. (5): Peak flux 
density of the emission. Col. (6): Deconvolved radius (at 5 cr contour level) of the cloud. Col. (7): Deconvolved FWHP line width. Col. (8): 
Integrated flux. Col. (9): Virial mass. Col. (10): Molecular gas mass from SMA C0(y=2-1). Col. (11): Surface density of the clouds, computed 
from the virial masses (Col. 9). Col. (12): Surface density of the clouds, computed from the molecular gas masses (Col. 10). 
" Deconvolved. 



To compute the expected millimeter continuum flux from the 
dust, associated with the CO clumps (with gas mass Mg), we can 
use the equation: 



(Launhardt & Henning 1997'), assuming optically thin thermal 
dust emission and isothermal conditions. By{Td) is the Planck 
function at the assumed dust temperature Tj, d the distance of 
the object and k{v) is the dust opacity per gram of gas. We adopt 
opacity equal to 0.005 cm^ g"' (Andre et al. 2000). For the de- 
rived masses of the individual CO clumps in Table [3] we esti- 
mate typical 1.3 mm continuum fluxes of 1-2 mJy (or 1-2 cr), for 
a range of temperatures between 20 K and 50 K. Therefore, it 
is not surprising that also from the CO clumps we do not detect 
1.3 mm continuum emission. 

Incidentally, we note that the continuum emission from dust 
may be a very good tracer of dense and massive molecular 
clouds. Indeed, an unresolved clump of mass Md may be more 
easily detected in continuum emission than in CO. This is be- 
cause, for a source smaller than the synthesised beam, one has 
for the ratio of continuum flux 5 d to CO flux S co that: 

S_A_ ^ K{v)Ma 

5 CO ttAvTJ^j 

where is the clump radius and Av is the CO line width. The 
latter in a virialised situation is proportional to (Mci//?ci)'^^ and 
so one concludes that the ratio of continuum to CO fluxes varies 
as R^^-^ for an unresolved clump. We assume in the above that 
the CO emission is optically thick and thermalised whereas the 
dust emission is optically thin. For a cloud of radius 10 pc, the 
ratio is almost ten times larger than the ratio for structures of 
the size of the beam {R^\ ~ 40 pc). This ratio can thus be sig- 
nificant if one is, for example, observing something similar to 
the "Massive Molecular Aggregates" postulated by Kobulnicky 

6 Johnson (2000), which are expected to have masses of or- 
der 10^ Mo and dimensions of order a few parsecs. Such objects 
would have millimeter fluxes of a few mJy and are barely out 
of reach for our observations. Clearly, verifying this speculation 



requires higher frequency, sensitivity and angular resolution sub- 
millimeter observations of the continuum from Henize 2-10. 

5. Discussion 

5.1. Dynamical properties of tlie moiecuiar gas 

An important aim of our study was to examine the kinematics of 
the molecular clouds for comparison both with the velocity field 
of the stars (see e.g. Marquart et al. 2007) and with earlier work 
on the gas kinematics by Kobulnicky et al. (1995). This is use- 
ful also for comparison with the molecular gas mass estimates 
derived in the previous section. 

Thus in order to investigate the rotation curve in the inner 
region of Henize 2-10 near the site of massive star formation 
and compare it with the results of Kobulnicky et al. (1995), we 
extracted a position-velocity diagram by taking a 4" wide slice 
of our SMA CO(/=2-l) data cube, passing through the brightest 
optical starburst region, at the same position angle as Kobulnicky 
et al. (1995) (PA=130°), which they find to be the direction of 
the steepest velocity gradient. The orientation of the position- 
velocity cut is shown in Fig|5] The position-velocity diagram, 
shown in Fig. |6l reveals two structures of which the innermost 
shows spatial velocity variations of 60 km s"' on a scale of 7 
arcseconds. This velocity gradient implies a dynamical mass of 
1.6 10^/(sinO^ Mq, within a radius of 150 pc, where / is the in- 
clination of the object with respect to the line of sight (/=90° 
means an edge-on rotating structure). 

Kobulnicky et al. (1995) detected in their OVRO CO(y=l-0) 
data a dynamically distinct feature, which was denoted by the 
letter "C". This CO feature, around 875 km s ' (see Fig. |6]l, is 
also detected in the SMA CO(7=2-l) integrated map between 
875 and 890 km s"', at the position 02001) = 8''36"'15!5 and 
^2000 = -26°24'39". Its SMA CO deconvolved size is about 
1'.'8 X 2'.'3, which corresponds to a linear size of 80 x 100 pc. Its 
total integrated flux density is 4 Jy km s"', consistent with the 
measurement by Kobulnicky et al. (1995), implying a molecular 
mass of 1.5 10* M©, using the same method described in sec- 
tion l4.2l (see Eq.[T]). Assuming a deconvolved velocity width of 
14.1 km s"' and a radius of 1", we estimate for this feature a 
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Fig. 6. SMA CO(/=2-l) position-velocity diagram of Henize 2- 
10 obtained from a 4" wide slice through the data cube, at a 
position angle of 130° (see Fig. |5]l, passing through the center 
of the map: 08'>36'" 15^.20, -26°24'34'.'0 (J2000). The horizontal 
axis represents the positional offset in arcsec with respect to this 
central position. Contour levels are from 25% to 95% in step of 
10% of the peak. 



virial mass of 1.8 lO^'M© (see Eq.|2]). This cloud is not listed in 
Table |3] since its peak flux density is below the detection thresh- 
old (15 cr, see ^4.1b that we have set for the output given by 
Clumpfind. 

It is interesting to compare the derived dynamical mass with 
the estimated stellar mass within the newly formed super star 
clusters and the mass of the molecular gas in the region. From 
our SMA CO(y=2-l), the mass of the molecular gas within this 
region is around 10^ Mq, which means that if the inclination of 
the molecular rotation axis is 90° (edge-on), the newly formed 
stars must contain at least 40% of the total mass in the inner 
150 pc radius. If the inclination is lower, the total dynamical 
mass is larger and, given that the total gasmass is determined by 
our SMA observations, the fraction of the total dynamical mass 
comprised by the stars will be larger. 

Conti & Vacca (1994) present Hubble Space Telescope im- 
ages which show 9 UV-bright knots in the inner 90 pc. They esti- 
mate a total mass of newly formed stars of 2.7 10^ M©. The upper 
limit to the stellar mass from our SMA observations is 6 10^ 
(see above), derived for a region which is about 3 times larger 
in radius than the one considered by Conti & Vacca (1994). This 
means that, extrapolating the estimated stellar mass from Conti 
& Vacca (1994) to the dimensions of our region, one obtains a 
stellar mass value between 2.7 lO^'M© and 2.4 10^ M©, respec- 
tively assuming that the stars are concentrated in the inner 90 pc, 
observed by Conti & Vacca ( 1 994), or uniformly distributed over 
a region which is 3 times larger in radius. Therefore, the picture 
presented in Conti & Vacca (1994) is consistent with the esti- 
mates which we discuss above. 

Another way to estimate the total mass involved in the whole 
region and how it is distributed between the molecular gas and 
the stars is to compute the (statistical) velocity dispersion, cry, of 
the clouds identified in the SMA CO(7=2-l) emission and derive 
the total mass of the region, assuming virial equilibrium. From 
Eq.|2l taking the FWHP velocity ( V8 /n(2) x cry), we derived a 
total mass of 1.2 10** Mq. Given that the total gas mass derived 
from our SMA CO(7=2-l) data is about 4-6 10^ Mq (derived 



from the molecular and virial masses respectively; see TableO, 
we estimate again that the stars should account for ~30-50% of 
the total mass, without taking into account the contribution from 
atomic gas, which is difficult to estimate at these angular scales. 

Our result is also in agreement with the OVRO CO(7=1-0) 
motions measured by Kobulnicky et al. (1995). Their study re- 
veals a rotation of the CO gas on large scales, with a velocity 
gradient which is about half of the gradient derived from our 
SMA CO data. They find that the low resolution of their obser- 
vations (6'.'5 X 5'.'5 beam size) gives only an "average" rotation 
curve smoothed over the beam and decreases the observed ve- 
locity gradient, which is consistent with what we observe. Their 
estimated dynamical mass is 3.2 10''/(sin/)^ Mq, which implies 
that the stellar component may contain up to half of the dynam- 
ical mass in the inner region, in agreement with our result. 

The gas velocity field is also consistent with Marquart et al. 
(2007). They find spatial velocity variations up to 60 km s"' for 
the ionized interstellar medium, but no systematic trend in the 
velocity of the stars and conclude that the stellar kinematics is 
governed by random motions. 

5.2. The physical properties of Clouds associated with SSCs 

Some views of cluster formation expect prestellar cores to have 
the same physical properties as the young stellar clusters they 
will eventually produce. Indeed, in the Galactic context the 
Infrared Dark Clouds (IRDCs), considered to be possible high- 
mass pre-stellar cores and revealed by absorption of the Galactic 
diffuse infrared background (Perault et al. 1996; Egan et al. 
1998), have similar surface densities as clusters associated with 
young massive stars in our Galaxy (Tan 2008). Therefore, in the 
extragalactic context, the expectation is to find extremely dense 
molecular clouds associated with SSCs. However, the limits im- 
posed by the resolution and sensitivity of the existing instru- 
ments make difficult the task of testing this scenario and resolv- 
ing the parent dense cores of the individual SSCs. A 5 cr detec- 
tion of an unresolved clump of 10^ Mq requires a sensitivity of 
~0.1 mJy beam"', assuming a temperature of 20 K, still a chal- 
lenge for the current millimeter interferometric facilities. 

The clouds, resolved in the CO(2-l) emission with the SMA, 
have masses and surface density values (see TableO of the same 
order of mag nitude of the galactic GMCs (M ~ W - 10^ Mq, 
2 ~ 0.035 g cm"^; Solomon et al. 1987), but systematically 
higher For a virialised cloud, surface density is a measure of 
pressure and it is comprehensible that the GMCs in Henize 2-10 
are at higher pressure than in the Galaxy. Moreover due to our 
limited angular resolution, our column density estimates may 
well be considerable underestimates. It seems reasonable there- 
fore that, as in the Milky Way, clusters in Henize 2-10 form from 
small parsec sized clumps at much higher pressure than the rest 
of the GMC. Using SMA, we have made a first step towards 
resolving the molecular clouds associated with the SSCs and to- 
wards testing for the presence of very dense and compact molec- 
ular clouds, predecessors of the SSCs. 

The density traced by the CO(7=2-l) line is limited (critical 
density n„ ~ 2.7 10^ cm"-') and to probe high density molecular 
gas (hht > lO"* cm"-'), which is expected to be associated with 
star (SSCs) forming regions, higher density tracers like HCO^ 
and HCN are needed. Indeed, Henize 2-10 is known to show 
strong HCO^ and HCN emission (Imanishi et al. 2007; our own 
IRAM 30 m spectrum presented in this paper), which reflect the 
spatial distribution of the dense molecular gas. The interferom- 
eter observations of Imanishi et al. (2007) have a beam size of 
5'.'5 X 10'.'8. To compare our SMA CO(7=2-l) with the NMA 
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Fig. 7. HCO+(7=1-0) integrated map from Imanishi et al. (2007) 
(beam size of 5'.'5 x lO'.'S, solid contours) and SMA CO(y=2-l) 
emission (dashed contours) smoothed to the resolution of the 
HCO^(7=1-0), overlaid on the VLA 3.6 cm continuum in grey 
scale (Johnson & Kobulnicky 2003). 



maps of Imanishi et al. (2007), we have smoothed our maps to 
the same angular resolution. The comparison is shown in Fig.|7] 
Given that the position calibration uncertainties for the interfero- 
metric data are much less than 1", the overlay of the CO(y=2-l) 
and HCO^(7=1-0) emissions seems to reveal a shift between the 
two molecular emissions, which is almost 2". This is surprising 
on one side because one would expect the HCO^(y=l-0), being 
a tracer of high density gas, to be associated to the peaks of the 
CO(y=2-l) molecular emission. On the other hand, the CO dis- 
tribution is complex in both position and velocity and it is likely 
that only some of the molecular clouds may be sufficiently dense 
to show strong HCO^(7=1-0) emission and this might produce 
the shift between the centers of the two integrated emissions. 
Besides, there might be also an effect due to the different opacity 
of the two lines. Finally, the HCO^(7=1-0) peak from Imanishi 
et al. (2007) appears to be closer to the SSCs in Henize 2-10 
than our SMA CO(7=2-l) peak. This could mean that the high 
density gas is mainly close to the SSCs. 

Indeed, the sensitivity and resolution of the HCO^(y=l-0) 
observations of Imanishi et al. (2007) are not adequate in space 
and velocity to resolve the dense gas peaks associated with the 
SSCs. It is thus necessary to observe high density tracers at 
higher resolution to resolve the molecular clouds out of which 
the SSCs formed. 



6. Summary and Conclusions 

In this paper we have presented new single dish and interfero- 
metric observations of the molecular gas and millimeter contin- 
uum of the starburst region in the dwarf galaxy Henize 2-10. 

Our observations confirm earlier detections of the '^CO and 
HCN(l-O) line emission. The relatively strong detections of high 
density molecular tracers associated with the young Super Star 
Clusters confirm that this galaxy is undergoing vigorous star 
formation and is an ideal laboratory to study extragalactic star- 
bursts. 



In this context, our CO(2- 1 ) and 1 .3mm continuum SMA ob- 
servations, with a linear resolution of 60 - 80 pc, represent a 
first attempt to resolve the parent molecular clouds out of which 
SSCs may form. We reveal a rich population of molecular clouds 
with estimated masses and densities in the upper range of those 
measured in our Galaxy. The molecular gas accounts for approx- 
imately half of the total mass in the inner region of the galaxy, 
while the young stellar clusters account for the remaining mass. 

We find possible evidence that the super star clusters are as- 
sociated with very massive and dense molecular cores, but our 
observations and tracers do not allow us to confirm unambigu- 
ously their presence. New higher angular resolution and sensi- 
tivity observations of high density tracers and (sub-)millimeter 
continuum are required to clarify this possibility. 
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